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CN ■ Abstracts. 

Qh' We have introduced and investigated so-called Schlomilchs and Bell's series for 

I modified Bessel's functions, namely, their asymptotic and non-asymptotic proper- 

ties, connection with Stirling's and Bell's numbers etc. 

We have obtained exact constants in the moment inequalities for sums of centered 
^ ■ independent random variables, improved their asymptotical properties, found lower 

^ . and upper bounds, calculated a more exact approximation, elaborated the numerical 

! algorithm for their calculation, studied the class of smoothing, etc. 
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O ■ 1 Introduction. Statement of problem. 

^ ' Let I,^{z) be usually modified Bessel's function: 

>< : W) = 2-^ z'- 4:~'z''/{k\ Tiu + k + 1)), (LO) 

where r(-) denotes the gamma-function and u > 0. 

We define the so-called Schlomilch's functions ^3,^2, Fi of the first kind for the 
values p>0, 6 > 0, (3 > : 

00 

k=—oo 



00 



F2{p;P)=Fs{p;l,P)= J2 l^l' W), (1.16) 



A:=— 00 



Fr(p) = F,{p, 1,1) = J2 \k\'h{l), (1.1c) 



A;=— 00 
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We also define the Schlomilch's functions of a second kind Gs,G2, G\ for positive 
integer values p— 1,2,3,... only and ^ > 0, /3 > as follows: 

oo 
fc=— oo 

oo 

G2ip;P)^Gsip,l,(3)^ Y: k^lM, (Lie) 

A;=— oo 
oo 

Gi{p)) = Gs{p, 1, 1) = Y (1.1/) 

A;=— oo 

Of course, G^{2m; 9, P) = F^{2m- 9,(3), m = 1, 2, 3, . . . . 
Recall that I-k{z) = Ik{z), k = 1,2, 3, ... Therefore, 

Fs{p;9,P)^Fs{p; 1/9, P). 

As for the definition in the case usual, i.e. non-modified Bessel's function and 
respective prehminary results see, for example, [1], [2], p. 38 - 42. 

Let us examine the probabilistic meaning of these functions. Let ^, r) he inde- 
pendent random variables characterized by the Poisson distribution with parameters 
A > and n > correspondently: Law{^) = Poisson{X), Law{rj) = Poisson[ii), 
or, in more detail: 

P(^ = A;) = exp(-A) X^/k\, A; = 0, 1, 2, ... ; 

P(,y = Z) = exp(-//) tx'/l\, Z = 0,1,2,...; 
and denote t = ^ — 77. Then we obtain for non-negative integer values n : 

00 

P(r = n) = 5^ P(^ = n + A;) P{r] = A;) = 

k=0 

00 

Y exp(-A - 11) n'' /{{k + n)\ n!) = 

k=i) 

exp(-(A + /x)) (A//x))"/2 /„ {2 , n = 0, 1, 2, . . . 

and therefore 

Elrr = exp(-(A + 11)) £ \nY (A//x)"/2 /„ (2 = 

n=—oo 

exp(-(A + /.)) F3 (p, ^/X/Ji, 2 , (1.2) 

00 

Er^ = exp(-(A + //)) ^ n^' (A//x)"/2 4(2 ^A^ = 

n=— 00 
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exp(-(A + //)) Ga (p, ^/x/Ji, 2 ^A^) . (1.3) 

(The case of negative integers n is considered in a similar way.) 
li X — II, then 

E|t|*' = exp(-2 A) F2(p,2 A), 
and in the case of X — /i — 1/2 : we obtain: 

E\t\p = Fi{p)/e. 

Here we introduce some new functions which represent generahzations of the 
classical Bell's number and functions. 

Generalized Bell functions of a first kind B^ij); a, X, 7)) can be defined as 

^oexp(A) •r(n + 7 + l) 
where A > 0, 7 > 0. We also define 

B,{p; a, A) = B,{p; a, A, 0) = ^ ' , (1.36) 

n=0 '<•• 

S2(p;A) = S3(p;0,A); (1.3c) 

S(p) = Brip) = S2(p; 1) = E (1-3^) 

We also define generalized Bell functions of a second kind L>4(p; a, A, 7)) for in- 
teger positive p values as follows: 

^ / ^ ^ ^ in-a)P X"" , , 

D,{p- a, A, 7 = E f\ J ^ ^-^^ 

^oexp(A) •r(n + 7 + l) 

A > 0,7 > 0. Wc also define Ds{p;a,X) = £)4(p; a, A, 0); D2{p; X) = 
D,ip;0,Xy, D,ip) = D2ip;l) = B^{p). 
The probabilistic interpretation is 

00 

E|{ - a\P = e-^ E 1^ - «rA7n! = B^ip; a, A) 

n=0 

and E(^-a)P = D3(l?;a,A). 

In this paper we investigate the estimations and asymptotic for intro- 
duced functions as p — > cxo. 

In Section 4 we find exact constants for moment estimations of sums 
of independent random variables. 
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2 Non-asymptotic properties. Relations with 
Stirling's numbers. 

Recall that the Stirling numbers of the second kind {s(n, r)} which have appeared in 
the combinatorial theory ([3], p. 117, [4], p. 234 - 243), are defined by the following 
identities: 

n 

= ^ s{n, r)x(^r)', X(r) — x{x — l){x — 2) . . .{x — r -\- 1), X(q) — 1. 

r=0 

1. Define the function Q2m{\ A*) for the values m = 1, 2, 3, . . . , A, // > by the 
following way: 



Q2m{\ =^ exp(-A - ii) Fs ^2m; \/x/Ji, 2\[\ii^ . 

Proposition: 

2m r 2m— r / 9m\ 

Q2m{X, /^) = EE E(-ir y 1^^ s{r,z)s{2m-r,j). 

r=Oj=0 j=0 \ ^ / 

Therefore, the function Q2m(AiA^) is a polynomial of the power 2m. 

Proof. As known Q2m{\ 1^) = — r/)^™". The conclusion 1 follows from the 
definition of Stirling's numbers, binomial formula and the following identity 

2. Let us consider the function D3(p;a,A) for integer non- negative values p. We 
have: 

Dsip; a. A) = E(e - af = E(-l)' f^) a^"^ E A^s(/, r) = 

1=0 \V r=0 

EA^E(-l)'f;V^-'«aO- 

r=0 l=r \V 

3. Let us denote for integer positive p the values: 

H(p; A, //) = exp(-A - //) G3(p; \/a7^, 2 ^A/I), 
where A = Ai + A2, A* = A*i + A*2, ^1^2 > 0, /xi,2 > 0. We assert that 

i7(p;Ai + A2,Aii + Ai2) = E T J^l'^; A^i) H{p - r; X2, fi2)- 

r=o vJ 

Namely, we can represent (on some Probability space) the random variables ^ 

and 1] as a sums C = + = ''7i + ^2, where all the variables are independent, 
Law{^i) — Poisson{Xi), Law{rji) — Poisson{iii), i = 1,2; then 
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H{p; A1+A2, //1+/X2) = E - 771) + (6 - V2)r = E f^) mi-viY m2-v2r-'' 

r=0 VJ 



r=0 



3 Asymptotic results. 

Denote q — q{p; a, 7, A) = (p — a — 7 — 0.5)/A, 

^ . ^ / X P . / X log log P 
elogp logp 

and assume a. 7, A = const, p 00. 

Let us introduce the following relations of equivalence: for two positive functions 
/2(p) we will write /i ~ /2 iff 

flip) ^f2{p){l + 0(1/ logp)). 

We can also write /i(-) x iff fl^'{p)/9iq) ~ /2/"(p)/^?(g). 

Note that if f{p) = C • h{p), C = const > 0, then /^^(p) ~ /iVp(p). 
Theorem 3.1. 

a. A-^sy^(p;a,A,7)Mg) ~ 1 + A(?). (3.1) 

b. Let us denote A = /5max(6', 1/9) for the function F = F2{p;9,(3). We assert 
that: 

F3(p;^,/3)xS4(p;0,A,0); (3.2) 

c.F,{p;9,P)>::Gs{p;9,py, (3.3) 

d.D4(p; a, A, 7) x B^ip; a, A, 7) ; (3.4) 
Proof. Step 1. We assert that A > =^ 

^<;,pVA)<^.f^!^. (3.5) 
n\ n\ A 



The first inequality is evident; let us prove the second one. We have: 

A^ 

fe=o^! (n + fc + 1)!" 



A-"4(2VA) = E 



1 A A^ 

. . < 



0!(n+l)! l!(n + 2)! 2!(n + 3)! 



Step 2. 



B^ip; a, A, 7) x A"53(p; A, a + 7). 
Proof. We can assume without loss of generality that a is integer. Further, 

, , ^ (k-ay X'' , ~ n^A" 
Let us denote for a > the function 

It follows from Stirhng's formula that at x > 1 

1 < < ^(«)- 



Hence 



A(q;) S4(p; a, A, 7) x = B^p - a - 7; 0, A, 0). 



m=l 

Thus, we have obtained a simple case i?2(-, ■) instead of B^l-). 

Step 3. We consider in this section the function F^{p]6,l3). We can suppose 
that ^ > 1. Denote 

00 

F+{p-e,(5)^Y.k''enM. 

k=0 

It follows, by virtue of the bide-sidc inequality < F3 < 2F~^ that F3 x F~^. 
Furthermore, we obtain using(3.5) 

F+ip; e,f3)^Y: -^rr^ ^ Mp; 0, a, O) = B,{p; A). 

k=l 

STEP 4. Instead of the general case, we need to investigate only the case B2. 
We obtain using the Stirling formula with remainder term again, B2 — -B2(p; A) = 

°° A"* °° 

X ^ exp (— mlogm + mlog(eA) + (p — 0.5)) . 

m=l m=2 

Replacing the last sum by an integral, we can see that 

B2 X exp(C/(p; x)) dx =^ V{p) = K 
U{p;x) — —X log X + X log(e A) + (p — 1/2) log x. 
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An integral for V can be estimated by means of classical saddle-point method 
[4], p. 119 - 127. 

In the next section we consider more exact estimations with non- asymptotic 
results in the particular case. 



4 Probabilistic applications and non-asymptotic 
results. 

Let p — const > 2, i = 1, 2, . . . , n be a sequence of independent centered 

E^(i) = random variables belonging to the space Lp, i.e. such that 

\\m\p = ^'^'mr<oo. m 

We denote Ea(i) = Er=ia(i), L{p) = CP{p), where 

C{p) = sup sup p \\^^i^)\\p (4_i) 

where "C" denotes the centered case, and the external " sup" is calculated over all 
the sequences of independent centered random variables that satisfy the condition 
(4.0). 

If all the variables in (4.0) are symmetrically distributed, independent, and 

||^(i)||p < oo, we denote the corresponding constants (more exactly, functions of p) 

def 

S{p) ("<S"'(-) denoting the symmetrical case) instead of C{p) and K{p) — S'^(p) 
instead of L{p). Obviously, S{p) < C{p), K{p) < L{p). It is proved in [5] that 
C{p) < 2S{p), L{p) < 2P K{p). 

The constant C(p), S{p) are called the exact constants in the moment inequalities 
for the sums of independent random variables. They play very important role in 
the classical probability theory ([6], 522 - 523, [7], p. 63;), in the probability theory 
on the Banach spaces ([8], [9], in the statistics and theory of Monte - Carlo method 
([10], section 5) etc. 

There are many publications on the behavior of the constants C{p),S{p) as 
p ^ oo. The first estimations where obtained in [11]; Rosenthal [12] proved, in fact, 
that C{p) < Cf ; Ci — const > 1; here and after Cj, j — 1,2,. . . are positive finite 
absolute constants. It is proved in [5], [13] that C{p) < 9.6 p/logp. In the works 
[14], [15] are obtained the non - asymptotic bide-sides estimations for S{p) : 

(eV2)-^ pI \ogp < S{p) < 7.35 p/ \ogp, p>2, (4.2) 

and there are some moment estimations for the sums independent non-negative 
random variables. See also Latala [5], Utev [16], [17]; Pinelis and Utev [17] etc. 

Ibragimov R. and Sharachmedov Sh. [8], [9] and Utev [17], [16] have obtained 
the explicit formula for S{p) : S{2) = 1; at p e (2,4] 



Sip) 
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P>4 ^5(p) = ||ti-T2|U (4.3) 

where a random variables Tj are independent and have the Poisson distribution with 

parameters equal to 0.5: Er,- = Var Tj = 1/2. 

As a consequence, it was obtained that as p ^ oo 

Si )= ^ (l \ (^ ^"^'^"^^ ^^ 
e-logpV V logP // 

The following representation for the values L{2m),m = 2, 3, 4, . . . : is obtained 
in [19]: 

oo 

C2'"(2m) = L(2m) = E(e - if^ = e"^ " l)^"*/^' = ^i(p), 

n=0 

where the random variable 9 has the Poisson distribution with parameter 1, and 
there is a hypothesis that for all values p > 4 C^(p) = L{p). 
Here we denote 



L(p)=E|^-l|f = e-^;^|n-l|7n! = B4(p; 1,1,0), G{p) = L^/P{p). (4.4) 

?i=0 

for all the values p > 4. In [20] estimations for C(p) are obtained in the case where 
the sequence is a sequence of martingale - differences, in [9], [8] are presented 

some generalizations for weakly dependent random variables {Hi)}- 

In this section we have improved the known bide-side estimations and 
asymptotic for 5'(p), G{p) at p — oo, found the exact boundaries for the 
different approximations of S'(p), G{p)\ describe the algorithm for the 
numerical calculation of K{p), L{p)] studied the analytical properties of 
K{p), L{p) etc. 

Note that there are many other statements of this problem: for the non-negative 
variables [12], [14]; for the Hilbcrt space valued variables ([8], [16]) etc. 
Let us denote at p > 4 : 5(p) = 1/ logp, 

Kp)=g{p) (l + A(p) + A2(p)) = 

[p/ (e logp)] ■ (l + log logp/ logp + (log logp/ logp)^) ; 
Theorem 4.1. 

1 = inf G'(p)/^(p) < sup G{p)/g{p) = C3, (4.5a) 

where 

Gs = supBl/^ip; l)/g{p) = G{G^)/g{G^) ^ 1.77638, 

p>4 
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C4 = argmax Bl^^{p)/g{p) ^ 33.4610; 

p>4 

(The equality ^ 1.77638 means that IC3 - 1.77638| < 5 • 10"^); 



1= inf Cip)/g{p)< sup C(p)/^(p) = Cs, (4.56) 

P=4,6,8... p=4,6,8... 



where C5 



where 



where 



inf l)/^(p) = G{C,)/g{C,) « 1.77637, Cg = 34; 

p=4,6,8,... 

1 = ini G{p)/h{p) < sup G'(p)//i(p) = G{C^)/h{Cs) = Cr, (4.5c) 

P>15 p>15 



C7 = sup 52''''(p; l)/h{p) « 1.2054, 

p>15 

Cg = argmaxS2/^(p; l)//i(p) 71.430; 

p>15 

1 = inf 5(p)/^(p) < s\x^S{p)/g{p) = Cg, (4.6a) 

P>4 p>4 



Cg = supiyi/^'(p)/5(p) = 5(Cio)/5(Cio) ^ 1.53572, 

p>4 

Cio = argmax Fl^^{p)/g{p) ^ 22.311; 
1 = inf 5(p)//i(p) < sup S{p)/h{p) = S{C,2)/h{Ci2) = Cn, (4.66) 

P>15 p>15 



where 



Cu = sup Fl'%p)/h{p) ^ 1.03734, 

p>15 

= argmaxFi/^(p)//i(p) ^ 138.149; 

p>15 

^ ^ ^M9n ^(Py^(P) < C{p)/h{p) = 

p=16,18,20,... p=16,18.20,... 

C(72)//i(72) = sup Fi/^(p;l)/%) fa 1.2053. (4.6c) 

p=16,18,20,... 

(We have choosen the value 15 as long as the function loglogp/logp monotoni- 
cally decreases at the values p > exp(e) ~ 15.15426). 
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Note that our estimations and constants (4.5a, 4.5b, 4.5c) and (4.6a, 4.6b,4.6c) 

are exact and improve the constants and estimations of Rosenthal [12]; Johnson, 
Schechtman, Zinn et al. [14], [15]; Ibragimov, Sharachmcdov [18], [19]; Latala [5]; 
Utev [16], [17] etc. For example, 1/(1/^2) ^ 1.41421, l.?>he/C^ « 11.2472. 

Theorem 4.2. At p — >• oo G{p) = \pj{e- logp)]x 



loglogp 1 log logp loglogp > 

1 + — I ^\ ^— ; — 2 ^-j — 2 — (l + o(l)) (4.7a) 

logp logp log> log> ' 

S{v) = [p/(e-logp)]x 



/ ^ log log P ^ 1 - log2 ^ log^ logp ^ ^ / log logp \ \ 
V logp logp log^p V log^p ]]' 

Let us denote hy N — N{p), M — M{p) the solutions of equations, which are 
unique: 

M{p) log M{p) = p, N{p) log(2iV(p)) = p, (4.8) 

for the values p > 4 such that N{p) = 0.5M{2p). 

Theorem 4.3. At p — >• oo, m = const = 2, 3, 4, ... , m — >• oo 

G{p) = M(p)^-^(f)/f exp(M(p)/p) (1 + 0(logp/p)), (4.9a) 



C(2m) =M^-^(2m)/(2m)(2m) exp(M(2m)/(2m)) (1 + 0(logm/m)), 

5(p) = N{p) {e/{2N{p)f^^/P (1 + 0(logp/p)). (4.96). 
Theorem 4.4. Let p be even: p = 2m, m = 2,3,4, ... . Then 

/r(2m) = , E E2-'-M2m-/,g).(/,r), (4.10a) 

1=0 \ * / q=0 r-=0 

C2"*(2m) = L(2m) = $^(-1)' s(2m -l,r). (4.106) 

i=0 \ W r=0 

For integer odds values p = 5, 7, 9, . . . we obtain the representation 

G'ip) = L{p) = (2/e) + ti-^f(t]^^(P-f')- (4-11) 



k=0 



Proof. First, we consider some auxiliary results. 

1. In the symmetrical case, for all the values p e [4, oo) we have: 



K{p) = (2/e) Y: 4(1) = Fi'\p). (4.12) 

n=l 
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Namely, we obtain from (4.3) for the values ti, T2 at n = 1, 2, . . . : 

oo o— /cr) — (n+fc) 

P(x.-r. = .) = e-'i:^^^ = /.(l)/e. 

2. On the basis of the equality (4.12) we can offer the numerical algorithm 
for K{p) investigation, calculation and estimation. To improve convergence rate of 
series (4.12) we can write: 

27r 4(1) = exp(cos(^)) cos(n0) dO, 

J —IT 

(see, for example, [21], p. 958, formula 5.) After the integration by parts, we obtain 
27r 4(1) = (-l)'^n-2'" r (exp(cos^))(2"^) cos(n^) dO, 

J — IT 

where m = 1, 2, Using the method of mathematical induction, we conclude: 

(expcos(^))(2™) =exp(cos(^)) P2m(cos(^)), 
where P2m{x) are polynomials of the degree 2m, which can calculated by recursion 

P2m+2{X) = (1 - X') (P// + 2P4(x) + P,^{X)) - 
X (pL{x) + P2m{x)) 

with the initial condition Po{x) — 1. Therefore, we obtain the following representa- 
tion for K[p) : 

oo 

TT e K{p) = n^'-^™ / expcos(^) P2m{0) cos{n9) d9. (4.13) 

n=l ''-'^ 

3. Corollary. For even numbers p — 2m, m = 1,2,3,... all the numbers 

K[p) = K{2m), L{p) = L{2m) are integer. 
In fact, it follows from Equation (4.12) that 

K{2m) = (tt e)-^ T g^^'^\e) cos(n^) dO^ 
n=i 

e-i(exp(cos^))(2™)(0) = (-1)"'P2^(1). 

It is easy to verify that all the coefficients of polynomials P2m{x) are integer; 
thus, the number P2to(1) is integer. 

The second conclusion of our corollary follows from Equation (4.10), as long as 
all the Stirling's numbers are integer. 

4. For example, K{<o) = 31, L(6) = 41. For non - integer values p we can use 
the method described above. We have obtained using a computer program: 
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TABLE 1 



p 


i^lPj 




P 




J-lpj 


2 


1 


1 


10.5 


14000.4 


A -I ("^ £^ 

41385.2 


4 


4 


4 


11 


30403.2 


98253.7 


4.5 


6.3358 


6.6712 


11.5 


67091.3 


236982 


5 


10.4118 


11.7358 


12 


150349 


580317 


5.5 


17.686 


21.538 


12.0 


341951.2 


1.44191E+006 


6 


31 


41 


13 


788891.0 


3.63328E+006 


6.5 


55.819 


80.5508 


13.5 


1.84518E+006 


9.27951E+006 


7 


103.22 


162.7358 


14 


4.37346E+006 


2.40112E+007 


7.5 


192.45 


337.176 


14.5 


1.04998E+007 


6.29176E+007 


8 


379 


715 


15 


2.55231E+007 


1.66888E+008 


8.5 


757.7 


1549.28 


15.5 


6.27927E+007 


4.47926E+008 


9 


1126.5 


3425.7358 


16 


1.56298E+008 


1.21607 E + 009 


9.5 


3015.0 


7721.29 


16.5 


3.93475E+008 


3.33839E+009 


10 


6556 


17722 


17 


1.00153E+009 


9.26407E+009 



TABLE 2 



p 


K(p) 


L(P) 


17.5 


2.57666E+009 


2.59791E+010 


18 


6.69849E+009 


7.36008E+010 


18.5 


1.75916E+010 


2.106E + Oil 


19 


4.66582E+010 


6.08476 + Oil 


19.5 


1.24952E+011 


1.77473E+012 


20 


3.37789E+011 


5.22427E+012 


20.5 


9.21603E+011 


1.55177E+013 


21 


2.53714E+012 


4.64999E+013 



5. Using a discrete analog of the saddle - point method ([4], p. 262 - 264), [18]), 
we find that 

M(p) = [p/logp] •(l + £(p)), 

where at p — > 00 

e{p) = A(p) + A2(p) - 8{p) A{p) (1 + 0(1)). (4.14) 

Hence 

7V(p) = [p/log(2p)]-(l+£(2p)) = 
12 



[p/logp]- [l + A(p) + A2(p)-5(p) A(p) (1 + log 2) (1 + 0(1)" . 
We define the following functions and constants for the values p > Pq — 700 

Ci4 = (1 - loglogPo/logPo)"' ~ 1.402365, 



C 
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;1 + 4A2(Po))^/^ + i1 0.928958, 



C(p) = log2/log(2p), 



£+(p) = A(p) + C,,A'{p), e_{p) = A(p) + Ci5A2(p), 

M+ = M+(p) = b/logp] • (1 + £+(p)), 

M_ = M_(p) = b/logp] • (1 + £-(p)), 

7V+(p) = b/(e • log(2p))] ■ (1 + £+(2p)), 

iV_(p) = [p/(e-log(2p))]-(l + £_(2p)). 
More exact calculation show us that for all the values p> Pq 

M4p) < M{p) < M+(p), N4p) < N(p) < N+(p). 

Namely, one can readly see that Wp > Pq ^ 

M_ log M_ < p = M log M < M+ log M+. 

6. Let us denote 

hix,p) = xP/r{x + 1), b2ix,p) = x7 (2^ r{x + 1)) ; 

V{x, p) — p log X — X log X -\- X, 

X{p) = V{M{p),p)/p = supV{x,p)/p, 

X>4: 

W{x,p) — plogx — xlogx + x{l — log 2), 
Y(p) = W(N(p),p)/p = supW(x,p)/p. 

x>4 

We have using the equality (4.14): X{p) = log[p/(e • logp)] + 



(4.15a) 



(4.156) 
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A(p) + Sip) + [log(l + eip)) - e{p) + A(p)£(p)]+ 

{5{p){e{p) - log(l + sip))} - 5{p)e{p) log(l + e{p))) = 

log[p/(e-logp)] +Xo(p), 
where at p > Pq ^2(p) < ^o(p) < ^i(p), ^i(p) 

A(p) + (^(p) + A{p)e+{p) + S{p)[e+{p) - log(l + e+{p))], (4.16a) 

A(p) + + [log(l + - £_(p)]- 

-S{p)s4p) log(l + £_(p)). (4.166) 
The function p — > Xi(p), p e [Pq, oo) is monotonically decreasing and 

exp(Xi(Po)) < 1.7563, Jiim Xi(p) = 0. (4.16c) 
In the same manner we obtain: Y{p) — log[p/(e • logp)] + Yq{p), where 

ro(p) = iog(i-C(p))-(i + £(2p))x 

[1 - C(p) - A(2p) + 5{2p) log(l + £(2p))] + 5(2p)(l + £(2p)) 
logg{p)+Yo{p), Y2{p) < Yoip) < Y^{p), 
Y^ip) A(2p) + 5{2p) + (1 + £+(2p)) • 5{p) log2/(l + 5{p) log2)+ 

£+(2p)[A(2p) + 5(2p)], (4.16d) 

Y2{p) =^ A(2p) + 6{2p) + £_(2p)[A(2p) + 6{2p)l (4.16e) 

where the function p — > yi(p), p e [Pi, oo). Pi = 10^ is monotonically decreasing 
and 

exp(ri(Pi)) < 1.442; Fi(p) | 0, p ^ oo; ^^^mY^ip) = 0. (4.17) 

7. Upper bound for L{p). We assume in this section that p > Pq = 700. 
Using the well - known Stirling's formula, we obtain for the values p > Pq, : 

e • L(p) — 1.5 = 6i(n — l,p) < / 6i (x, p) (ia; + sup 6i(a;,p) < 

n=3 x>3 
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(27r)-^/^ exp(p • X{p)) + {27r)-^/^ exp{V{x,p)) dx. 
Splitting the last integral into three parts so that 



CO 



J{p) =^ exp(V(a;,p)) dx = Ji + J2 + J3, Jj = Jj{p): j = 1, 2, 3, 
where 

Ji{p) = / exp(y(x,p)) dx, J2 = / cxp{V{x,p)) dx, 
J2 Jm-^ 

roo 

Js^ ex.p{V{x,p)) dx, 

we obtain for the integral J2, taking into account the inequalities M_ < M < M4 
and inequality: x e [M — ^fp, M + ^] ^ 

V{x,p) < pX{p) - 0.5 (x - Mf ■ < 
pX{p) - 0.5(x -Mf -p- M-^{p) : 



the following: 

-VP 



h < exp(p • X{p)) ■ / exp (-0.5p{x - Mf Ml^) dx < 

exp(p • X{p)) ■ r exp (-0.5p(x - M{p)f M"^) = 
2^ exp(p • X{ip)) M+/^ < exp(p • X{p)) ■ 



where 



= -^2^ • [1 + A + CuA^] / logp. 



Now we estimate the integral J3. For the values x > M -\- ^ the following 
inequalities are valid: 

V{x,p) < pX{p) - 0.5 • (2p) • (p/MKp)) < pX{p)- 
log2p-(l + A + Ci4A2)'2. 

dV{x,p)/dx < -p/MKp) [x - M{p) - ^] ; 
therefore J3 < exp(p ■ X{p)) ■ ^'2(p), where ^'2(p) = 

exp (- log^p • (1 + A(p) + Ci4A2(p))-2) X 



15 



/ cxp ( -p M~^(x - M - a/2p) ) dx^ 
exp (- log^p . (1 + A(p) + Ci4A2(p))-2) X 

p-(l + A + Ci4A2)'-log-'(p) 

Analogously, we find the upper estimate for Ji. 
Thus, L{p) < ■ exp{p ■ X{p))x 

1.5exp(-p • X{p)) + (27r)-^/2 ^ ^^^^^ ^ 2(27r)-'/' ^2(p) 



exp(p-X(p)) ■ ^f(p), (4.18a) 
where we find by direct calculations: \I'3(Po) < 1.00826 and at p > Pq 

Mp) i 1, P ^ oo; *3(p) < 1 + Ci8 logp/p. (4.186) 

9. Lower bound for L{p). Denote q — p — 1/2. using Sonin's estimate for 
factorials, we obtain: 

oo oo 

eL{p) > J2 bi{n -l,p)^Yl ^Mp) > 



n=4 



n=3 



/•oo /•oo 

6i(x,p) > (27r)-^/^ exp(-l/12) e^p{V{x,q)) dx > 
(27r)-^/2 exp(-l/12) exp(y(x, ?)) dx. 



Since the following implication holds: x e [M(g), M(g) + y/g] 
y(x, q) > q X{q) - 0.5{x - M{q)f qMZ\q). 



we have: 



eL{p) > (27r)-^/2 exp(-l/12) exp{qX{q))x 



M(q) 



exp 



-0.5q MZ^{q) {x - M{q)f dx > M^{q) x 



0.5 exp(-l/12) ^ exp{qX{q)) [l - exp (-^7^-) 



where 



e-exp(p-X(p)) •*^(p). 



*4(p) i l,p ^ oo; *4(p) > 1 + Ciglogp/p. 
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Thus, 

exp(p • X{p)) ■ < L{p) < exp(p • X{p)) ■ (4.19a) 



*3(p) < 1 + Ci9 logp/p, *4(p) > 1 + C20 logp/p, (4.196) 

^3(p) i l,p ^ 00; ^3(^0) < 1.00826. (4.19c) 

10. Upper and lower bounds for K{p) are obtained analogously to the 
upper bound for L{p), but we assume in this section that p> Pi = 10^. In brief, 

~ 4-^ 1 ~ A-^ _ ^ 1.285 

^^k\ {n + k)\ ^ n! ~ ^ ^ n\ ' 

hence 

K{p) < 2e-^/^ X 2^ • ^^(p; 1/2). 

n=l ^• 

Further, we conclude, using the Stirling estimate for factorials again, that: 
0.5 e^'^ K[p) = '^b2{n,p) < / b2{x,p) dx + supb2{x,p) < 

n=l ^>2 

(27r)-^/2exp(p- F(p)) + {2tt)-^/'^ exp(iy(a;,p)) dx. 
Again, we spht the last integral: 

h j exp(iy(a;,p)) dx = 

( rN(p)-^ fN(p)+y/p roo \ 

/ + / + / exp(Vr(a;,p))o?x = 

75 + /g + h- As long as at X > N{p) + ^ =^ 

W{x,p) < pY{p) - 0.5p^Nl(p) = pY(p) - 0.51og'p • (1 + £+(2p))-^ 

dW/dx < -pN-\x - - Vp), 

we obtain: 

I7 < exp(py(p)) -p log-'p (1 + £+(2p))2x 

exp(-0.51og'p(l + £+(2p))-2). 
Further, if x e [iV(p) - y^, N{p) + y^] then 

1^(0;, p) < pY{p) - 0.5pAr-2(p) . _ ]sj(^p^y_ 

17 



Therefore, 

h < exp(py(p)) • (1 + £+(2p))/log(2p) 
and K{p) < 2e-^/^ exp(py(p))x 

(27r)-^/' + VP(1 + £+(2p))/log(2p) + 2(27r)-^/2p log-'p' 



X 



;i + e+i2p)Y ■ exp(-0.5 log'p (1 + e+{2p))-^ 
Lower bound for K{p). We obtain: 0.5 e K{p) > 



(4.20) 



5^n^2-'^/n! = 52(p;l/2) > exp(-l/12) (271)-^/^ exp(gF(g)) 



X 



n=l 



JNiq) 



exp 



-0.5 g Nz\q) {x - N{q)y 



dx > 



exp(-l/12) ^ q-'/' N_{q) (l - exp (-gV7V!(g))) / log(2g) = 



exp(-l/12) ^7r/2 exp(gr(g)) (l + £-(2g))x 

(l_exp {-qyN'_{q)))/\og{2q). 

Further estimations are similar to the estimation of L(p) and can be omitted. 
As a result, 



exp(p . Y{p)) ■ ^l(p) < K{p) < exp(p • Y{p)) ■ *f (p). 



(4.21a) 



where at p > Pi 



*5(p) < 1 + C21 logp/p, ^'6(p) > 1 + C22 logp/p; 



(4.216) 



*5(p) i 1, p ^ 00; *5(i'i) < 1.000833. 



(4.21c) 



11. For exact computations, we need to estimate the derivatives of our functions 

L(p), K{p). We show here the estimation of derivatives L^'") (p), m = 1,2 Namely, 

e-L("^)(p) = 



E 

n=3 



(n - ly log"*(n - 1) ^ g (n - 1)^ log"*n ^ 



00 



n 



n=2 



n\ \n>3 n 



We estimate the derivative K'^'^\p)^ m = 1, 2, ... in an analogous way. 



(4.22) 
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It follows from this estimation that the functions L{p) and K{p) are infinitely 
different iable at the interval p E (4, oo). As long as L(4 - 0) = L(4 + 0) = 4, K{4: - 
0) = K{4: + 0) =4, both functions K{-), L(-) are continuous in the semiclosed 
interval [2, oo). However 

^(4 - 0) = ^(4 -0)^ 3.149195, ^(4 + 0) ^ 3.51934, 
dp dp dp 

^(4 + 0) ^ 3.86841, 
dp 

and therefore, both the functions K{-), L(-) are not continuously differentiable 
in the set (2, oo). In the open intervals (2,4) and (4, oo) all the functions 
L{p),K{p),C{p),S{p) are infinitely differentiable (see (22) and [18], [19], [16] ). 

5 Proof of the probabilistic results. 

Proof of theorem 4.1. We find by direct calculations that G{Ci) / g{Ci) a; 1.77638, 
but for we conclude from (4.17) that for the values p> Pq — 700 

G{p)/g{p) < 1.00826 • 1.75913 = 1.77366, 

hence 

aigm.a.'xG{p)/ g{p) G [4,700]. 
pe[4,oo) 

We obtain by direct calculations using known numerical methods and by means of 
computer: 



Further, 



^max^^G(p)/^(p) = G{C,)/g{C,) « 1.77638. 



inf G'(p)/^(p) = min |^ min^^ G{p)/g{p), \rdG{p)/g{p) 



Our computations show that 

min G'(p)/^(p) 1.332, 

pe [4,700] ' 

and it follows from (4.18a), (4.18b), (4.18c) and (4.19) that 



Thus, 



^^l^G{p)/g{p) = jim G(p)/^(p) = 1. 



MG{p)/g{p) = lim^G{p)/g{p) = 1. 



Analogously, S{Cio)/g{Cio) 1.53572, but for the values of p > Pi the following 
holds: 
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S{p)/g{p) < 1.0008333 • 1.443 < 1.4444. 

Therefore 

argmax5(p)/^(p) e [4,1000000]. 

pG[4,oo) 

We have copmuted the following: 



Further, 



max = ,S(Cio)/y(Cio) ^ 1.53572. 

pe[4,Pi] 



inf = inin{ min inf S{p)/g{p)} 

p>4 Pe[4,Pi] p>Pi 



\YmS{p)/g{p)^l. 

Other assertions of theorem 1 are obtained analogously. 

Proof of theorems 4.2 and 4.3. It follows from inequalities (4.19a), (4.19b), 
(4.19c) and (4.21a), (4.21b), (4.21c) that 



exp(X(p)) • (1 + C20 logp/p) < G{p) < exp(X(p)) • (1 + C19 logp/p), (4.23a) 



exp(r(p)) • (1 + C22 logp/p) < Sip) < exp(r(p)) • (1 + C21 logp/p). (4.236) 



Substituting the expression (4.16a) and (4.16b) into equation (4.23a), we obtain, 
after simple calculation, our assertions (4.9b). Wc obtain (4.9a) in a similar way. 

Finally substituting expressions (4.16d,c) into (4.23b), we obtain (4.7a), (4.7b). 

Proof of Theorem 4.4 Since E6'(^) = 1, r = 0, 1, 2, . . . , we conclude 

k k 

E9'' = E^s(A;,0^(r) = ^^^(fc,/). 
1=0 1=0 

Equation (4.11a) follows from the binomial formula. The equality (4.10) is proved 
analogously. 

Let us prove (4.11b). Since 



^n7n! = e-5i(p)-l, p=l,2,3,...; ^ 1/n! = e - 2, 

n=2 n=2 

we conclude for the values p = 5, 7, 9, . . . : e • L(p) = 

n=2 "'■ n=2 k=0 V^/ 



1 + 



E(-l)'Q(e5i(p-A;)-l) 
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(eS,(0) - 2) = 2 + e • ±{-1^ f^) B,{p - k) - £(-1)'= (( 

k=o V^/ fe=o V^y 

2 + eE(-l)'f>Vi(p-^)- 



k=0 



6 Some generalization and concluding remark. 

Our results allow us to obtain some generalizations of symmetrically distributed ran- 
dom variables {r]{i)}, i = 1, 2, 3, on the Hilbert space. Let {H, 1 1 1 • 1 1 1) be a sep- 
arable Hilbert space with the norm ||| • |||, P(?7(0 & H) — 1, yi — 1,2,3 . . . H^yllp 
EVf(|||r^(i)|||f)<oo, p>4, 

IIE^WIIp 



Z(p) — sup sup ■ 



n"max(||Er?WI|2,(E|ki)||^)V^)- 
Utev ([16], [17]) has proved that Z{p) = S{p), p > 4 (in our notations). Therefore 

1 = inf Z(p)/^(p) < supZ(p)/^(p) = Cg ^ 1.53572, 

P>4 p>4 

1 = inf Z{p)/h{p) < sup Z{p)/h{p) = Cu ^ 1.03734. 

P>15 p>15 

Moreover, let us examine the following problem. For the values A,De (0, oo), p — 
const > 4 we define 



QP{p, A, D) = sup sup E|||V^( 
n>l {e{i)} i=i 



^)lll^, 



where the interior " sup " is calculated over all the sequences of if — valued in- 
dependent symmetrically distributed random variables {0{i)} under the following 
conditions: 

n n 
i=l i=l 

We denote 

i = 0.5(^/L'P)V(P-2), 
then t e (0, 1/2]. Utev [16], [17] has also proved that 

QP{P, A, D)^(A- D~^y^^~'^ E\u^ - v^\^, 

where v\, vi are independent Poisson-distributed random variables with parameter 
t : 

P(u. = k)= exp(-t) • t^/k\. A; = 0, 1, 2, ... . 

Hence 
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n=l 



We also denote for p > 4, t e (0, 1/2] 

R{p,t) = Q{p,A,D)-(DyA) 



l/(p(p-2)) 



2e-^* £n^4(2t), 



n=l 



It = u{t) — sup R{p,t)/g{p), T — T{t) — argmax F{p,t), 

p>4 p>4 

and using the above-described method, that at p — > oo, t = const e (0, 1/2] 
R{p, t) = [iM(p/t)]i-*^(^/*)/^ (^1 + o (^i^ j j , 



loglogp 1 + logt log^logp /loglogp' 
1 H ; ^ ; \ ^ — 9 \- o 



logp 



logp 



log^p 



log p 



i=0 \ W q=0 r=0 

mmp,t)/g{p) = lim i?(p,i)/^(p) = 1; 

p>4 p— >oo 



s{l,r) 



TABLE 3 



t 


T(t) 


u(t) 


0.45 


26.228 


1.48566 


0.4 


32.206 


1.43438 


0.35 


42.120 


1.3815 


0.3 


60.67 


1.3281 


0.25 


102.47 


1.2732 


0.2 


145.96 


1.2163 



Note that at t = 1/2 => T(l/2) = do ^ 22.311, m(1/2) = 1.53572. 

Apparently, it is interest to obtain exact constants in the moment inequalities 
for sums of independent non-negative random variables in the spirit of Refs. [12], 
[14], [5] etc. 
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